I. Introduction

Energy is the cornerstone of modern economies, not only for the suppliers of physical goods and services, but also as a means of social welfare and comfort for people in general. Thus, an understanding of how price changes impact the energy demand of suppliers and consumers is very important (Gil-Alana et al., 2020). In the last year, energy prices in countries such as Spain and Portugal rose from approximately €50/ MWh in November 2020 to approximately €200/MWh by the end of October 2021, causing inflation tensions along the different elements of the supply chain. In European countries, the energy pricing strategy is fixed with day-ahead market systems (single-day-ahead coupling).

In the Spanish market, aligned with the objectives of the European Internal Energy Market (Directive 2009/72/EC), the Operador del Mercado Ibérico de Energía (OMIE) is the nominated electricity market operator for managing the day-ahead and intraday electricity markets by submitting selling and takeover bids with market-based supply prices. This market has a large number of regulatory specificities and competition constraints compared to other European energy markets. The main issues to be resolved in the Spanish market are the reduced competition with few suppliers, the tariff deficit due to past regulated prices below production cost, potential market integration between different players, and the introduction of new technologies to reduce CO2 emissions (Duarte et al., 2017).

Additionally, from a global perspective, a final challenge involves the interconnection networks with foreign countries as a protective measure. Currently, the European Union (EU) energy policies called for a well-integrated internal energy market in 2020 by achieving the interconnection of at least 10% of the installed electricity production capacity for all EU member states, with a 15% target in 2030 (Rubino & Cuomo, 2015). However, in the specific case of the Spanish market, these interconnections can be done only across the Pyrenees mountains, and thus there is a real risk of Spain become an energy island. An adequate cross-border interconnection capacity should help to avoid the internal development of dispatchable reserve capacity, for balance and grid security purposes. N. C. Figueiredo et al. (2015) tracked the Spanish case, showing that Spain has already surpassed this value, reaching 25.6%, and is aiming to achieve 3,000 MW in the near future, which will represent 32% of the maximum demand considered in this study; however, by 2020, overall consumption was still locally fulfilled.

Several authors have recently investigated the Spanish market and its peculiarities. In terms of the green renewal process R. Figueiredo et al. (2019) studied the replacement of traditional coal-based power plants with photovoltaics, while Pereira and Saraiva (2013) explained the implications of the penetration of wind power. Ciarreta and Zarraga (2010) studied the dynamic relation between electric consumption and the gross domestic product, and, more recently, Ciarreta and Zarraga (2016) examined the volatility of the intra-day electricity market. Rubino and Como (2015) explained interconnection needs, Capitán Herraiz and Rodriguez Monroy (2013) studied the futures market, Gil-Alana et al. (2020) examined the relation between energy consumption and prices in both the futures and spot markets, and Lagarto et al. (2014) studied the market power of generating firms in the day-ahead market.

The objective of this paper is to evaluate the impact of shocks in the Spanish energy market by studying the degree of persistence in hourly energy prices during the last 12 months, to understand if shocks that are affecting this market are expected to remain permanent. From a methodological perspective, fractional integration will be used as a time series technique, allowing for a more flexible dynamic specification of the data time series than other approaches based on integer degrees of differentiation.

II. Data description and methodology

The data of this paper are based on 12 independent time series associated with the hourly pricing for each of the last 12 months (November 2020 to October 2021), with hourly records (00:00 to 23:00) of the spot day-ahead market price traded by OMIE. These series include both the shock of Storm Filomena last winter and the recent natural gas rally shocks. These 12 series of data provided by OMIE are plotted in Figure 1 and their descriptive statistics summarized in Table 1.

Figure 1
Figure 1.Detail of day-ahead energy pricing of Spanish Market over the last 12 months with the different time series under analysis.

Notes: Timeseries are built with hourly samples for each day for the specific month. For example, October 2021 has 24 hourly samples x 31 days = 744 samples. Source: Authors’ using data from OMIE.

Table 1.Descriptive statistics of the datasheet
Series Mean Price
(EUR/MWh)
Std. Dev.
(EUR/MWh)
Maximum
(EUR/MWh)
Minimum
(EUR/MWh)
November 2020 41.91 9.05 62.38 8.00
December 2020 41.97 12.57 68.90 1.95
January 2021 60.17 26.89 121.24 0.16
February 2021 28.49 19.46 80.00 0.16
March 2021 45.45 15.67 85.43 4.04
April 2021 65.02 15.35 93.92 0.90
May 2021 67.12 21.17 97.01 0.01
June 2021 83.30 12.55 105.51 3.84
July 2021 92.42 12.62 111.36 2.67
August 2021 105.94 17.63 137.46 0.90
September 2021 156.14 24.50 221.00 99.99
October 2021 199.90 48.58 319.03 15.00

This table shows descriptive statistics of the series. Std. Dev. denotes standard deviation.

As far as the methodology is concerned, we use fractional integration, which is widely used in the analysis of energy prices. The applied model is the following:

y(t)=α+βt+x(t),(1B)dx(t)=u(t),t=1,2,...

where y(t) refers to the observed data; α and β are unknown parameters referring to an intercept and a linear time trend; B is the backshift operator, that is, Bx(t)=x(t1), and x(t) is integrated of unknown order that is estimated from the data; and u(t) is an uncorrelated zero-mean process.

This model is applied separately for each monthly datasheet, yielding 12 independent results. Here, the estimation of the differencing parameter d is crucial. Thus, if d=0, x(t)=u(t) in (1), and x(t) is said to indicate short memory, as opposed to the case of long memory that takes place when d>0. From a statistical viewpoint, the cutoff point is 0.5. Thus, if d<0.5, x(t) is covariance stationary; however, it becomes nonstationary for d0.5 and even more nonstationary as we increase the value of d, where we note that the variance of the partial sum increases in magnitude with d. Finally, from a policy perspective, mean reversion occurs if d<1, and shocks will have permanent effects if d1.

III. Empirical results

The results of this analysis appear to be quite similar over the 12 monthly series, with d values ranging from 1.30 to 1.53, providing empirical evidence that all the series under study are highly persistent. Table 2 displays the estimates of d in equation (1) under the assumption of white noise errors, capturing all time dependence with the differencing parameter d. Thus, there is clear evidence of non-mean reversion, since d is greater than one in all cases, as well as in all confidence intervals. Thus, shocks affecting these series are expected to have permanent effects.

Table 2.Estimates of d in the model given by equation (1)
Estimates of d under different scenarios
Series No terms* An intercept An intercept and a time trend
November 2020 1.38 (1.27, 1.50) 1.38 (1.27, 1.49) 1.38 (1.27, 1.50)
December 2020 1.38 (1.29, 1.48) 1.53 (1.43, 1.64) 1.53 (1.43, 1.64)
January 2021 1.46 (1.37, 1.57) 1.53 (1.42, 1.64) 1.53 (1.42, 1.64)
February 2021 1.48 (1.35, 1.61) 1.48 (1.35, 1.62) 1.48 (1.35, 1.62)
March 2021 1.30 (1.18, 1.43) 1.36 (1.24, 1.49) 1.36 (1.24, 1.49)
April 2021 1.26 (1.11, 1.38) 1.43 (1.30, 1.56) 1.43 (1.30, 1.56)
May 2021 1.35 (1.24, 1.46) 1.47 (1.35, 1.59) 1.47 (1.35, 1.59)
June 2021 1.17 (1.08, 1.28) 1.39 (1.26, 1.51) 1.39 (1.26, 1.51)
July 2021 1.17 (1.08, 1.29) 1.33 (1.20, 1.48) 1.33 (1.20, 1.48)
August 2021 1.26 (1.14, 1.37) 1.40 (1.28, 1.53) 1.40 (1.28, 1.53)
September 2021 1.16 (1.06, 1.27) 1.30 (1.18, 1.42) 1.30 (1.18, 1.42)
October 2021 1.29 (1.17, 1.41) 1.44 (1.31, 1.58) 1.44 (1.31, 1.58)
Estimated coefficients on the selected models above
Series D (95% band) Constant (t-value) Time trend (t-value)
November 2020 1.38 (1.27, 1.49) 19.285 (7.61) ---
December 2020 1.53 (1.43, 1.64) 47.865 (16.20) ---
January 2021 1.53 (1.42, 1.64) 52.207 (10.08) ---
February 2021 1.48 (1.35, 1.62) 2.329 (14.21) ---
March 2021 1.36 (1.24, 1.49) 36.356 (7.74) ---
April 2021 1.43 (1.30, 1.56) 56.596 (13.99) ---
May 2021 1.47 (1.35, 1.59) 71.098 (13.57) ---
June 2021 1.39 (1.26, 1.51) 89.969 (22.48) ---
July 2021 1.33 (1.20, 1.48) 94.156 (18.91) ---
August 2021 1.40 (1.28, 1.53) 77.329 (15.45) ---
September 2021 1.30 (1.18, 1.42) 132.569 (15.76) ---
October 2021 1.44 (1.31, 1.58) 224.606 (15.58) ---

*: " No terms" means that no deterministic terms are included in the model. The values in parenthesis in the upper part refer to the 95% confidence band for the values of the differencing parameter d, and those in bold correspond to the selected specification according to the deterministic terms. In the lower part, in parenthesis (in column 3) the t-values for the estimated constant.

However, these results are in contradiction with the findings of Gil-Alana et al. (2020), who analyzed future pricing markets and found d values around 0.62 with confidence intervals between 0.58 and 0.68. Given that our study is conducted with daily observations, but at longer periods of analysis with higher frequencies, further research should clarify if this d is modified by using longer periods of analysis (longer periods of data) or by increasing the sampling rates (hourly or daily).

IV. Conclusions

This paper has examined the hourly structure of 12 months of energy price data in Spain by using a long memory model based on fractional integration. The results indicate that all series are very persistent, with values of the differencing parameter substantially and significantly higher than one, implying a lack of mean reversion and the permanency of shocks. These results, however, could be due to the high frequency of the data used in the application. Robustness tests were conducted allowing for autocorrelation in the error term, and the results were very similar to the main results of the paper, finding support for high levels of persistence. Nevertheless, further research should be conducted on this issue. Finally, because the data include the period of the COVID-19 pandemic, the high degree of persistence observed in the data indicates that the present shock could have a permanent effect on the series, implying the need for strong action to recover the original long-term projections.


Acknowledgment

Luis A. Gil-Alana and Miguel Martín-Valmayor gratefully acknowledge financial support from grant PID2020-113691RB-I00 funded by MCIN/AEI/10.13039/501100011033. Comments from the editor and an anonymous reviewer are gratefully acknowledged.